Holes in a Quantum Spin Liquid 
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pLn ■ Magnetic neutron scattering provides evidence for nucleation of antiferromagnetic droplets around 
impurities in a doped nickel-oxide based quantum magnet. The undoped parent compound contains 
a spin liquid with a cooperative singlet ground state and a gap in the magnetic excitation spectrum. 
Calcium doping creates excitations below the gap with an incommensurate structure factor. We show 
that weakly interacting antiferromagnetic droplets with a central phase shift of vr and a size controlled 
by the correlation length of the quantum liquid can account for the data. The experiment provides 
I a first quantitative impression of the magnetic polarization cloud associated with holes in a doped 
^ transition metal oxide. 

Spin density modulations in transition metal oxides are receiving huge attention because of possible connections 
to high temperature superconductivity. The modulations appear upon introduction of charge carriers, via chemical 
substitution, into an insulating and antiferromagnetic parent compound, and tend to be static when the carriers are 
frozen, dynamic when they are mobile. Evidence for such modulations has been largely confined to materials whose 
magnetism and charge transport are quasi-two-dimcnsional(i-^), and whose parent insulators are ordered antiferro- 
magnets. We provide evidence for analogous phenomena in a quasi-one-dimensional oxide(5), Y2_j;Ca2,BaNi05, for 
Q I which the parent is a spin liquid by virtue of quantum fluctuations JO), and which has motivated considerable 
Q ■ theoretical activity (i J- J 6). 

The key features of this orthorhombic material are the chains of NiOg octahedra (Fig. lA). The octahedra are 
' corner-sharing, which results in the dominance of magnetism((5) and electrical conduction (5) by the very simple 
^ ' ..0-Ni-O-Ni-O.. back-bone. The magnetic degree of freedom at each Ni site is the spin S=l associated with the 3d^ 
ly-j ■ configuration of Ni^+. Each of the S—1 ions is coupled to its neighbors via antiferromagnetic (AFM) super- exchange 
' through shared ions. Replacing the off-chain Y^^^ by Ca^+ introduces holes primarily onto apical oxygen atoms in 
the chains and while the materials remain insulators, doping significantly increases conductivity at finite temperatures. 
Thus, Y2-iCa:rBaNi05 is a one-dimensional analog of the cuprates, where off-Cu02-plane chemical impurities donate 
holes to the Cu02 planes. 

Magnetic one-dimensionality causes the parent compound Y2BaNi05 to be a spin liquid prevented from ordering 
antiferromagnetically by quantum fluctuations. The material is not an ordinary paramagnet with heavily damped 
' spin fluctuations, but rather the magnetic analog of super-fluid ^He because it has a macroscopically coherent quan- 
tum ground state and lacks static order. Affleck, Kennedy, Lieb and Tasaki (AKLT) provided a prescription for 
this sid.te{17,18). They considered each S—1 degree of freedom as the triplet ground state formed between two fer- 
romagnetically coupled S=l/2 degrees of freedom, in accord with atomic physics which posits a large and positive 
Ph , Hunds rule coupling between the two S=l/2 holes of Ni^+. To understand the qualitative features of the ground and 
^ ' excited states of the spin-1 chain, they started from the limit of vanishing intra-atomic interactions between spin- 1/2 
holes. Because each site of the spin chain has as many nearest neighbors as S=l/2 degrees of freedom, there exists 
a macroscopic singlet with the full symmetry of the Hamiltonian wherein all S=l/2 degrees of freedom take part in 
nearest neighbor bond singlets (Fig 2A). Excitations from this "Valence Bond Solid" (VBS) state are triplets formed 
by broken valence bonds that propagate along the chain. The triplet band is separated from the VBS by an energy 
of order J corresponding to the energy required to create a triplet at rest. As the Ferromagnetic (FM) Hund's rule 
coupling is taken to much larger values than J these low energy states evolve adiabatically into the ground and excited 
states respectively of the S=l chain. The major change is that for energies below the intra-atomic exchange energy, 
spin-1/2 degrees of freedom can be observed only as edge states, not in the bulk. 

The triplet band can be observed experimentally for Y2BaNi05 (Fig. 2A). An image of the magnetic neutron 
scattering intensity as a function of wave vector transfer q along the spin chains (in units of the inverse lattice 
spacing) and energy transfer, Tiuo, was obtained using time-of-flight neutron spectroscopy at the ISIS facility in the 
UK. Magnetic neutron scattering occurs only when huj and q lie on a dispersion relation, thus providing evidence for 
coherent triplet propagation along spin chains. What distinguishes this excitation from a spin wave is that it propagates 
in a medium with no static magnetic order. Key features are (i) a 9 meV Haldane gap(_Z5), which corresponds to the 
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threshold energy for triplet creation, (ii) the magnetic excitation bandwidth of 62 meV that translates into an AFM 
exchange constant J = 21 meV, and (iii) vanishing intensity for 17= n27r, which is evidence for a singlet ground state. 

What happens when by Ca doping we add carriers to the quantum spin liquid of Fig. 2A?. Our x = 0.095(5) 
sample (Fig. 2B) reveals excitations following a similar dispersion relation as for pure Y2BaNi05. However, near the 
former minimum at = tt, their intensity is substantially reduced and new magnetic scattering appears below the 
Haldane gap. More detailed data collected near and below the gap (Fig. 3B) show that near the gap enc;rgy, the 
scattering is clearly peaked at g = tt but for fiw = 3 — 7 mcV it is most intense along two vertical lines, displaced 
symmetrically about tt. Higher statistics energy-integrated (/-scans are shown (Fig. 4) that accurately quantify the 
incommensurate peaks for three single crystals with different levels of doping. Surprisingly, increased doping docs not 
yield a proportional increase in the shift of the peaks from tt. Specifically, if we fit the three data sets by the sum of 
two identical Lorentzians centered aX -k ±5q with half- widths-at- half-maximum, k, we find that 8q extrapolates to a 
finite value in the dilute limit increasing with x as Sq{x)/'iT « 0.059(2) + 0.25(3)a; (sec inset to Fig. 4(B)). A similar 
x-dependence is found for k{x)/Tr « 0.047(3) + 0.20(6)x with the implication that a single length scale, which is finite 
in the limit a; — »■ 0, controls both q{x) and k{x). 

There are various models which could produce the incommensurate fluctuations we have discovered in 
Y2-i;Caj:BaNi05. The flrst proceeds from a charge-ordering hypothesis where the holes order periodically to minimize 
Coulomb repulsion and the spin system acquires a period commensurate with that of the hole lattice. The magnetic 
fluctuations are spin waves broadened in q due to disorder in the hole lattice. Such a scenario is commonly supposed 
for the two-dimensional analog, doped La2Ni04, where key experimental support came from electron(^O) and neutron 
diffraction(^i,i3i?). Thus inspired, we used both techniques to search for hole ordering in Y2-2;Ca2;BaNi05 down to 
15K. The failure of the search makes the charge ordering explanation of our experiments unlikely. Furthermore, in 
the simplest picture, the hole density has the same period as the spin density squared so that 5q{x)/'K « x, which is 
entirely inconsistent with Fig. 4 (dashed line in the inset). 

A second possibility is to identify the incommensurate magnetic fluctuations with electron-hole pair excitations 
from the one-dimensional hole liquid associated with mobile S=l/2 hole propagation through the VBS(iJ). The 
wave-vectors it ±5q would be vectors spanning the Fermi surface while the vertical nature of the incommensurate 
streaks (Fig. 3) would imply a Fermi velocity exceeding 0.5 eV-Angstrom. The problem with this picture is that even 
if we can ignore that the samples seem not to be proper metals - the resistivities (5) are high at low T - the Fermi sea 
is expected to grow in proportion to x leading to 5q{x) « x, a result again contradicted by Fig. 4. 

A last possibility starts from consideration(i5) of the magnetic structure factor for a single static hole donated to 
the NiO chain by the Ca ions. Because it carries a spin, a hole on a super-exchange mediating oxygen atom induces 
an effective ferromagnetic Ni-Ni interaction(^5,.2.^). The hole-carrying oxygen divides the chain into two semi-inflnite 
segments, for which wave- functions can be constructed from a basis set which is the direct product of wave- functions 
for each segment and for the hole on the intervening oxygen. The ground state wave-function for a semi-infinite 
segment is a doublet rather than the singlet wave-function for the unterminated chain. Pictorially, (Fig. IB-C) the 
chain ends have unpaired spins because only one of the two AKLT spins, which together represent the S=l sites, can 
be paired with an AKLT spin on a neighboring site. The effective ferromagnetic interaction between semi-infinite 
chains yields a ground state with inversion symmetry about the oxygen site. The edge states are not localized on Ni^+ 
atoms neighboring the hole, but extend into the chain with AFM modulation over a distance of order the Haldane 
spin correlation length, As a result, ground state correlations {SoSj) between spins at site and site j, where 

the FM bond occurs between sites —1 and 0, should take the form (— l)^e~'-'l"^ for j > and (— l)J+ie~l-'+^l'* for 
j < -1- 

Neutron scattering measures the modulus squared Fourier transform of the spin polarization. For the ground state 
associated with a single FM bond inserted within an infinite chain, this is 

(l + e'-)cos(g/2) ^ ^ 
cosh K + cos q 

Eq. 1 describes a function with nodes for q = (2n -I- l)7r and pairs of peaks, with half widths and incommensurability 
of order k (see dashed line in Fig. 4A). Thus, the description leading to Eq. 1 docs have a chance of accounting for a 
key feature of our data. There are however three difficulties, the first being that the nodes predicted by Eq. 1 are not 
seen experimentally. The second problem is that we have neglected the spin on the intervening oxygen. The third is 
that Eq. 1 describes the ground state for a single FM bond, implying that neutron scattering should be elastic, not 
inelastic as seen experimentally. 

The first difficulty has a simple resolution. At finite hole densities the polarization clouds overlap and the isolated 
impurity model is inadequate. A crude model, which considers the overlaps, simply truncates the polarization clouds 
at neighboring impurity sites. Because FM impurity bonds are randomly distributed, the inversion symmetry char- 
acterizing isolated impurities is broken thus allowing intensity at g = tt. Because we know the impurity density, x, 



S{q) = 



2 



from neutron activation analysis, the only parameters in such a description arc the extent of the polarization cloud, 
K~^, which we adjusted to optimize the fit to our data. As shown by the red lines in Fig. 4 the model provides a good 
account of the data with k^^ = 8.f (2), 7.3(2), and 7.2(5) for x = 0.04,0.095, and 0.f4 respectively. These values are 
close to the exponential decay length of 6.03 calculated for the AFM spin polarization at the end of an S=l chain(^5). 

The modeling described so far does not include the spins of the holes responsible for the effective FM couplings 
between Ni^+ ions. The holes reside in oxygen orbitals of Y2_2:Caj,BaNi05(5), but are almost certainly not confined 
to single, isolated oxygens. We have consequently generalized Eq. 1 to take into account the hole spins, with - for 
the sake of definiteness - the same net amplitude as either of the Ni^+ spins next to the FM bond, and distributed 
(with exponential decay) over £ lattice sites centered on the FM bond. As long as £ exceeds the modest value of 2, 
comparable to the localization length deduced from transport data (5), the pronounced asymmetry about tt which 
occurs when ^ = is relieved sufficiently to produce fits indistinguishable from those in Fig. 4. 

How do we account for the inelasticity of the incommensurate signal? One approach is to view the chain as consisting 
not of the original S—1 degrees of freedom but of the composite spin degrees of freedom induced around holes. The 
latter interact via overlapping AFM polarization clouds and hole wave functions, to produce effective couplings of 
random sign because the impurity spacing can be even or odd multiples of the Ni-Ni separation. With weak interchain 
coupling the ground state is likely to be a spin glass, as deduced from other experiments (S6') on Y2_xCaxBaNi05. 
The "incommensurate" nature of the excitations continues to follow from the structure factor of the spin part of the 
hole wavefunctions. 

In summary, we have measured the magnetic fluctuations in single crystals of a doped one-dimensional spin liquid. 
At energies above the spin gap, the triplet excitations of the parent compound, Y2BaNi05, persist with doping. How- 
ever, below the gap, we find new excitations with a broad spectrum and characteristic wave vectors that are displaced 
from the zone boundary by an amount of order the inverse correlation-length for the parent. The incommensurate 
fluctuations, encountered here in a doped one-dimensional transition metal oxide, resemble those seen in metallic 
cuprates. However, one-dimensionality makes them easier to model for Y2BaNi05, and our analysis reveals that "in- 
commensurate" peaks arise naturally even without hole order because of the characteristic spin density modulation 
that develops around a defect in the singlet ground state of a quantum spin liquid. This phenomenon accounts for the 
weak dependence of the incommensurability on doping in the one-dimensional nickelate. Indeed, Y2_2:Ca2;BaNi05 
gives us our first quantitative impression of the magnetic polarization cloud associated with the holes in a doped 
transition metal oxide. Our results imply that the spin part of the hole wave-function is actually the edge state 
nucleated by the hole in a quantum spin fluid. 
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FIG. 1. (A) Chain unit of Y2BaNi05 featuring nickel atoms with octahedral oxygen coordination and Ca^+ impurities on 
sites. (B) Schematic of S= 1/2 de grees of freedom at the edges of the AKLT state. (C) Ferromagnetically coupled chain-end 
S=l/2 degrees of freedom. 

FIG. 2. Overview of magnetic fluctuations in (A) pure and (B) 9.5% doped Y2BaNiOr, at T=10K. The initial neutron beam 
energy was 90 meV and the chain axis was perpendicular to the incident beam direction. Boxes indicate regions examined in 
Fig. 3. The color bar shows values for audEf units of mbarn meV~^ per Ni. 

FIG. 3. Low energy detail of magnetic excitations in (A) pure and (B) 9.5% doped YgBaNiOs. (A) shows time-of-flight 
data at T=10K (MARI spectrometer, ISIS pulsed neutron source) while (B) shows data collected at T=1.5K using triple-axis 
spectrometers (SPINS and BT2 at NIST steady state nuclear reactor with final energies 5 meV and 14.7 meV respectively). 
The color bar shows values for ^ dcidE ™ units of mbarn meV~^ per Ni. 



FIG. 4. g-scans, collected using SPINS at NIST, through the incommensurate peaks at fixed energy transfer of 4.5 meV for 
three different hole concentrations. The energy resolution of the spectrometer was 2 meV FWHM. The dashed green line in A 
shows the single impurity model Eq. 1 convolved with the instrumental resolution (solid bar). The red lines take into account 
that neighboring impurities truncate the spin polarization cloud around an impurity bond. Inset in B shows half the distance 
5q{x) /tt between the peaks of two Lorentzians superposed to fit the data. 
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